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Abstract: In this work we obtain the DLLA limit of BFKL in the 
dipole picture and compare it with HERA data. We demonstrate that in 
leading - logarithmic - approximation, where a s is fixed, a transition between 
the BFKL dynamics and the DLLA limit can be obtained in the region of 
Q 2 w 150 GeV 2 . We compare this result with the DLLA predictions obtained 
with a s running. In this transition is obtained at low Q 2 (< 5 GeV 2 ). 

This demonstrates the importance of the next-to-leading order corrections 
to the BFKL dynamics. Our conclusion is that the F% structure function is 
not the best observable for the determination of the dynamics, since there is 
great freedom in the choice of the parameters used in both BFKL and DLLA 
predictions. 
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1 Introduction 



The behaviour of ep/pp scattering in high energy limit and fixed momentum 
transfer is one of the outstanding open questions in the theory of the strong 
interactions. In the late 1970s, Lipatov and collaborators |JJ established the 
papers which form the core of our knowledge of Regge limit (high energies 
limit) of Quantum Chromodynamics (QCD). The physical effect that they 
describe is often referred to as the QCD pomeron, or BFKL pomeron. The 
simplest process where the BFKL pomeron applies is the very high energy 
scattering between two heavy quark-antiquark states, i. e. the onium-onium 
scattering. For a suficiently heavy onium state, high energy scattering is a 
perturbative process since the onium radius gives the essential scale at which 
the running coupling is evaluated. Recently |^|, |[ || this process was studied 
in the QCD dipole picture, where the heavy quark-antiquark pair and the 
soft gluons in the large N c limit are viewed as a collection of color dipoles. 
In this case, the cross section can be understood as a product of the number 
of dipoles in one onium state times the number of dipoles in the other onium 
state times the basic cross section for dipole-dipole scattering due to two- 
gluon exchange. In Mueller demonstrated that the QCD dipole picture 
reproduces the BFKL physics. In this work we discuss the BFKL pomeron 
using the QCD dipole picture. 

Experimental studies of the BFKL pomeron are at present carried out 
mainly in HERA ep collider in deeply inelastic scattering in the region of low 
values of the Bjorken variable x = ^— , where Q 2 = — q 2 . Here p is the four- 
momentum of the proton and q is the four-momentum transfer of the eletron 
probe. In this case, QCD pomeron effects are expected to give rise to a 
power-law growth of the structure functions as x goes to zero. However, this 
study in deeply inelastic scattering is made difficult by the fact that the low x 
behaviour is influenced by both short distance and long distance physics [[J . 
As a result, predictions at photon virtuality Q depend on a nonperturbative 
input at a scale Qq < Q. This makes it difficult to desantangle perturbative 
BFKL predictions from nonperturbatives effects. Moreover, the program of 
calculating the next-to-leading corrections to the BFKL equation was only 
formulated recently |J. Of course there are uncertainties due to subleading 
corrections and from the treatment of infrared region of the BFKL equation 
which will modify the predictions of this approach. 

One of the most striking discoveries at HERA is the steep rise of the 
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proton structure function F2(x,Q 2 ) with decreasing Bjorken x [[[]. The be- 
haviour of structure function at small x is driven by the gluon through the 
process g — > qq. The behaviour of the gluon distribution at small x is itself 
predicted from perturbative QCD via BFKL equation. This predicts a char- 
acteristic x~ x singular behaviour in the small x regime, where for fixed a s 
the BFKL exponent A = ^-Aln2. It is this increase in the gluon distribu- 
tion with decreasing x that produces the corresponding rise of the structure 
function. However, the determination of the valid dynamics in the small x 
region is an open question, since the conventional DGLAP approach || can 
give an excellent description of F2 at small x. The goal of this letter is to 
discuss this question. 

Recently Navelet et al. [[J applied the QCD dipole model to deep inelastic 
lepton-nucleon scattering. They assumed that the virtual photon at high 
Q 2 can be described by an onium. For the target proton, they made an 
assumption that it can be approximated by a collection of onia with an 
average onium radius to be determined from the data. This model described 
reasonably the F 2 data in a large range of Q 2 (< 150 GeV 2 ) and x. 

In this letter we obtain the double-leading-logarithmic-approximation 
(DLLA) limit of BFKL in the QCD dipole picture using the approach pro- 
posed by Navelet et al.. This limit is common to BFKL and DGLAP dy- 
namics. We show that using our DLLA result HERA data can be described 
in the range Q 2 > 150 GeV 2 and all interval of x. Moreover, we compare 
our results with the predictions of Navelet et al. and with the predictions 
obtained using DGLAP evolution equations in the small-x limit. This letter 
is organized as follows. In Section 2 the QCD dipole picture is presented. In 
Section 3 we obtain the proton structure function in this approach and its 
DLLA limit. In section 4 we apply our result to F2 HERA data and present 
our conclusions. 

2 QCD dipole picture 

In this section we describe the basic ideas of the QCD dipole picture in 
the onium-onium scattering. Let A be the scattering amplitude normalized 
according to 
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The scattering amplitude is given by 

A=-if d 2 x 1 d 2 x 2 J dz 1 dz 2 $ ( -°\x l , z x ) $ (0) (x 2 , z 2 )T{x_ x ,x 2 ) , (2) 

where zi) is the squared wave function of the quark-antiquark part of 

the onium wavefunction, x^ being the transverse size of the quark-antiquark 
pair and Z{ the longitudinal momentum fraction of the antiquark. In lowest 
order T is the elementary dipole-dipole cross-section odd- 

In the large N c limit and in the leading-logarithmic-approximation the ra- 
diative corrections are generated by emission of gluons with strongly ordered 
longitudinal momenta fractions z, t » Zi+\. The onium wave function with 
n soft gluons can be calculated using perturbative QCD. In the Coulomb 
gauge the soft radiation can be viewed as a cascade of colour dipoles ema- 
nating from the initial qq dipole, since each gluon acts like a quark-antiquark 
pair. Following [0, [|, we define the dipole density n(Y,x,r) such that 

N(Y,r) = [ dzi I d 2 x$ {0) (x,z 1 )n(Y,x } r) (3) 



is the number of dipoles of transverse size r with the smallest light-cone 
momentum in the pair greater than or equal to e~ Y p + , where p + is the light- 
cone momentum of the onium. The whole dipole cascade can be constructed 
from a repeated action of a kernel K, on the initial density n (x,r) through 
the dipole evolution equation 

n(Y, x, r) = n (x, r) + / dy dsJC(r,s)n(y,x,s) . (4) 

JO JO 

The evolution kernel K, is calculated in perturbative QCD. For fixed a s 
and in the limit iV c — > oo the kernel has the same spectrum as the BFKL 
kernel. Consequently, the two approaches lead to the same phenomenological 
results for inclusive observables. The solution of (^) is given by [| 

1 xexp[(a P - 1)Y] ( nln 2 (x/r) \ 
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where ap — 1 = (8aCV/7r)/n2. 

The onium-onium scattering amplitude in the leading-logarithmic approx- 
imation will be written as in (|2]), but where T is now given by 

/(I T cP S 
-=-=n(y/2,^ 1) r)n(y/2,x 2 ,i)a DC . (6) 

Consequently, the cross section grows rapidly with the energy because the 
number of dipoles in the light cone wave function grows rapidly with the 
energy. This result is valid in the kinematical region where Y is not very 
large. At large Y the cross section breaks down due to the diffusion to large 
distances, determined by the last exponential factor in (|5|), and due to the 
unitarity constraint. Therefore new corrections should become important 



and modify the BFKL behaviour 1C 



The result ([5]) was obtained for a process with only one scale, the onium 
radius. In processes where two scales are present, for example the eletron- 
proton deep inelastic scattering, this result is affected by non-perturbative 
contributions |J. Therefore, the application of BFKL approach at ep scat- 
tering must be made with caution. 



3 Structure function in the QCD dipole pic- 
ture 

Our goal in this section is to obtain the proton structure functions 
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using the QCD dipole picture. In order to do so we must make the assumption 
that the proton can be approximately described by onium configurations. 
Basically, we make use of the assumption 

afp = atr™ x V , (8) 



where V is the probability of finding an onium in the proton. In order to 
obtain the <Jtot° mum we wm follow IJ, where the factorization [11 
used in the context of the QCD dipole model. We have that 



was 



a 



7 onium 
L.T 



J d 2 rdz^ 0) (r,z)a^ dipole (x,Q 2 J r) , (9) 
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where the 7* — dipole cross section reads 

gV*^ e (x,Q 2 ,r) = J d 2 k £^^9(^)8(z,k,d • (10) 

In (0) ay g is the 7*g — > qq Born cross section and Q is the non-integrated 
gluon distribution function. The relation between this function and the 
dipole density is expressed by 

/s 2 f ' dzt z 
— \ — n ( z, iL,s)a rd ( — ,s 2 k 2 ) , (11) 
s Jo zf zf 

where n(z/, r, s) is the density of dipoles of transverse size s with the smallest 
light-cone momentum in the pair equal to z/p + in a dipole of transverse size 
r, of total momentum p + . This is given by the solution (||). 
After some considerations Navelet et al. obtain 

n J Lixi 7 

(12) 

where ^^xil) is BFKL spectral function (for more details see 0). 

The 7* onium cross section is obtained using flT2"| ) in (|D|). The result 
depends on the squared wave function <3>(°) of the onium state, which cannot 
be computed perturbatively. Consequently, an assumption must be made. 
In H this dependence is eliminated by averaging over the wave function of 
transverse size 

J dzd 2 r(r 2 y&°\r,z) = (M 2 ) -7 , (13) 

where M 2 is a scale which is assumed to be perturbative. Therefore the 
proton structure functions reads 

FlH*,?) - ^ J ^y^Vy^v-Vv^tfi) , (14) 

where P(7,M 2 ) is the Mellin-transformed probability of finding an onium 
of transverse mass M 2 in the proton. Using an adequate choice for this 
probability (see ref. 0) the proton structure function can be written as 

fi,x(*,Q 2 ) = ^ /^x,r(7)^(gr^*>'»^( 7 ) • (15) 

7T J Lixi 7 Un 
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The expression fll^D can be evaluated using the steepest descent method. 
The saddle point is given by 



X'ds) = • (16) 

TT X 

Using in the expression ([16]) the expansion of the BFKL kernel near 7 = |, 
we get 



Consequently 



^7C(3)/ni 



F 2 (x,Q 2 ) = F T (x,Q 2 ) + F L ( Xl Q< 



where 



1 Q [ ( 4aAr nM ); l_a; n 2 Q] 

Ca 2 —e iy n x 2 Q o' , (18) 
Qo 



(19) 



7C(3)/r»i 



The parameters C, Q and a are determined by the fit. Using C = 0.077, 
Qo = 0.627 GeV and a ~ 0.11, Navelet et a/, obtained that the expression 
flT|) fits the HERA data in the region Q 2 < 150 GeU 2 . 

In this letter we analize the behaviour of F2 obtained by expression (|15|) in 
the double leading logarithmic approximation (DLLA). This limit is common 
to both DGLAP and BFKL dynamics, i.e. it represents the transition region 
between the dynamics. Therefore, before the determination of the region 
where the BFKL dynamics (BFKL Pomeron) is valid, we must determine 
clearly the region where the DLLA limit is valid. In this limit x{l) = ■ 
Using this limit in (|16D we get that the saddle point is at 

Is = • (20) 
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Consequently, we get 



The result (|21|) reproduces the behaviour of double leading logarithmic ap- 
proximation. As this result was obtained considering the dipole model, then 
a s is fixed. The parameters C, a and Qq should be taken from the fit. In the 
next section we compare this result with HERA data. 



4 Results and Conclusions 

In this section we compare the expression (]21f) with the recent HI data [^]. 
In order to test the accuracy of the F 2 parameterization obtained in formula 
fl2lD, a fit of HI data has been performed. The parameters obtained were 

C = 0.0035 , Q = 0.45 and a = 0.19 . (22) 

In the figures ([I]) and (|2|) we present our results at different Q 2 and x. 
The predictions of Navelet et al. (dashed curve) are also presented. While 
the expression (|TiJ) describes the data in region Q 2 < 150 GeV 2 , we can 
see that the expression ([H]) (solid curve) describes HI data in kinematical 
region Q 2 > 150 GeV 2 . Therefore the HERA data are described by the 
DLLA expression in the high Q 2 region. Moreover, we can conclude that 
there is a transition between the BFKL behaviour and the DLLA behaviour 
in the region Q 2 « 150 GeV 2 . This could be the first evidence of the BFKL 
behaviour in F2. However, this conclusion is not strong since that (|18|) and 
( |21~D were obtained in the leading-logarithmic-approximation, where a s is 
fixed. Moreover, the parametrizations obtained using the DGLAP evolution 
equation describes the HERA data |T2[. 

The program of calculating the next-to-leading corrections to the BFKL 
equation is not still concluded []nj. However, some results may be antici- 



pated. For instance, the NLO BFKL equation must have as limit the DLLA 
limit in the region where a s log -log Q 2 ~ 1. This limit is common to 
both DGLAP and BFKL dynamics. As the NLO corrections to the DGLAP 
evolution equations are known, the DLLA limit with a s running is well estab- 
lished. The DLLA limit obtained considering the DGLAP evolution equation 
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was largely discussed by Ball and Forte Consequently we can estimate 
the importance of a s running in our result. In figure (Q) we compare our 
results with the DLLA predictions with a s running (dot-dashed curve), ob- 
tained using DGLAP evolution equations in the small-x limit. In this case 
the DLLA limit can describe one more large kinematical region. The region 
Q 2 < 5 GeV 2 is not described by DLLA a s running. Consequently, the kine- 
matical region where the BFKL dynamics may be present is restricted to the 
low Q 2 region. 

Our results are strongly dependent on free parameters, since there is 
great freedom in the choice of parameters used in both BFKL and DLLA 
predictions. This comes from theoretical uncertainties, for example, the de- 
termination where the pQCD is valid (i.e. the Qo value). However, our 



qualitative result agrees with the conclusion obtained by Mueller [14], that 



demonstrated that the BFKL diffusion leads to the breakdown of the OPE 
at small-x. Using the bounds obtained by Mueller we expect that the BFKL 
evolution can be visible in a limited range of x at Q 2 — > Qq. Moreover, our 
result agrees with the conclusion of Ayala et al. ||15|| , where the shadowing 
corrections to F-i structure function were considered. In this case the anoma- 
lous dimension is modified and the BFKL behaviour only can be visible in 
the region of low Q 2 . 

In this paper we calculate the DLLA limit of BFKL in the dipole pic- 
ture. The determination of this limit is very important, since it is com- 
mon to BFKL and DGLAP dynamics. Therefore, before the determina- 
tion of the region where the BFKL dynamics is valid, we must determinate 
clearly the region where the DLLA limit is valid. We demonstrate that in 
the leading-logarithmic approximation (a s fixed) a transition region between 
the BFKL dynamics and the DLLA limit can be obtained in the region of 
Q 2 w 150 GeV 2 . We compare this result with the DLLA predictions ob- 
tained with a s running. In this transition region is obtained at low Q 
(< 5 GeV 2 ). This demonstrates the importance of the next-to-leading order 
corrections to the BFKL dynamics. Our conclusion is that the F 2 structure 
function is not the best observable in the determination of the dynamics. 
From the inclusive measurements of F2 it seems improbable to draw any 
conclusion based on the presently available data. It is theoretically question- 
able whether it will be possible as long as one considers only one observable. 
The better observables for determination of the dynamics are the ones asso- 
ciated with processes where only one scale is present, since in these processes 
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it is inambigously possible to isolate the effects of the BFKL behaviour. 
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Figure Captions 



Fig. |I|: Behavior of proton structure function predicted by BFKL (|18|) 
(dashed curve) and DLLA (|21| ) (solid curve). Data of HI JJ]. See text. 

Fig. 0: Behavior of proton structure function predicted by BFKL (|TE|) 
(dashed curve), DLLA (|2l| ) (solid curve) and DLLA with running coupling 
constant (dot-dashed curve). Data of HI 0]. See text. 
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Figure 2: 
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